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Abstract. In this paper we compare two polymer stretching experiments. The 
outcome of both experiments is a force-extension relation. We use a one- 
dimensional model to show that in general the two quantities are not equal. In 
certain limits, however, both force-extension relations coincide. 



1. Introduction 

Traditional experiments on biological molecules are carried out on large groups 
of molecules. Such groups of molecules are macroscopic thermodynamic systems 
which can be interpreted using conventional thermodynamics. Recently, it became 
possible to experimentally manipulate individual molecules. Such single molecules are 
investigated for several reasons, some of which are mentioned below. 

The mechanical properties of single molecules are important in processes like 
polymer folding 0Ej and DNA strand separation 

Biophysicists study small motors that are responsible for converting chemical 
energy to useful forms of work. These biological machines work under non-equilibrium 
conditions. The working principles are poorly understood. Many biomolecules can 
function as such motors. So the study of single molecules can help to understand the 
properties of biological machines EJ ■ 

The experiments [HI G] of interest for the present paper combine results obtained 
under varying thermodynamic circumstances. In one type of experiment, the two 
ends of a single linear molecule are held at fixed positions and the exerted forces are 
measured. In the other type of experiment, the force exerted on the ends is constant 
and the position of the end points is measured. The corresponding ensembles differ 
because in one case the end positions are fluctuating quantities while they are kept 
constant in the other case. Nevertheless, the outcomes of the two experiments should 
be related, as argued in =6. . In general, results of statistical physics do not depend on 
the choice of ensemble. This equivalence of ensembles has been proven for large classes 
of systems. However, it may fail in case of single molecule experiments for at least two 
reasons. The proof of ensemble equivalence relies on the thermodynamic limit, while 
here the number of degrees of freedom is limited. It may also break down because 
long range interactions, in particular those due to excluded volume, are important. 
The present paper is based on the model studied in This model is simple enough 
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so that statements can be made which do not depend on the assumption of a large 
number of degrees of freedom. 

In the next section we discuss the ensemble dependence of measurements on single 
molecules in more detail. In section[3]we give a short review of the results obtained in 
jS] and we show how this results can be applied on single molecules. In sections 0] and 
El we derive expressions for the force-extension relation for two different ensembles, 
the canonical and the grand canonical ensemble. In section we compare the force- 
extension relations of both experiments. The final section gives a short discussion of 
the results. 

2. Single molecule experiments 

In the next two paragraphs we discuss two types of experiment, first from experimental 
and then from theoretical point of view. 

The first experiment is the ideal isotensional experiment where the forces on the 
end points of the molecule are kept constant and the fluctuating end-to-end distance 
is measured. The second experiment is the ideal isometric experiment where the end 
points of the molecule are held fixed and the fluctuating force is measured. A clear 
description of a possible experimental setup can be found in 9 . The outcome of both 
experiments is a force-extension curve. In the isotensional experiment, one obtains 
the average end-to-end distance (x)(F) as a function of the applied force F. In the 
isometric experiment, one obtains the average force (F)(x) as a function of the fixed 
end-to-end distance x. Inverting (x)(F) results in F((x)). An interesting question is 
whether (F)(x) and F((x)) should coincide. 

From a thermodynamic point of view, the end-to-end distance x is a (non- 
conserved) observable and the force F is the corresponding control parameter (just 
like temperature is the control parameter of the energy). By measurement of the 
average value of the observables one can estimate the value of the corresponding 
control parameters. In this view, the canonical ensemble is a one parameter family. 
Such family has only one fluctuating observable (mostly the energy) and one control 
parameter (mostly temperature). Hence, the isometric experiment is described by a 
canonical ensemble because only the energy fluctuates. The isotensional experiment 
requires a grand canonical ensemble because both the energy and the end-to-end 
distance fluctuate. One can derive expressions for (F)(x) and F({x)) in both 
ensembles. Whether the expressions are equal or not, boils down to a question of 
equivalence of a canonical and a grand canonical ensemble. 

In [H] Keller et al investigate the ensemble dependence of single molecule 
measurements. They replace the familiar thermodynamic potentials (Gibbs or 
Hclmholtz free energy) by appropriate potentials of mean force. Keller et al find that 
the two experiments should be related by a Laplace transform. In the thermodynamic 
limit (infinite long chain) the Laplace transform is dominated by it's saddle point value 
and the two experiments are then related by a Legendre transform [7] . Following this 
reasoning, the force-extension relations of both experiments coincide. However, in 
general the saddle point approximation is not valid. Therefore the experimentally 
obtained force-extension relations are distinct. 
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Figure 1. Two walks with 6 steps which end in x = 0, have the first step to 
the right and have 4 kinks, e and /i are the probabilities to go straight on when 
walking to the right, respectively to the left, po is the probability that the first 
step goes to the right. 



3. Model 

In [S] we have studied a non-Markovian random walk in one dimension. It depends 
on two parameters e and /x, the probabilities to go straight on when walking to the 
right, respectively to the left. We have calculated the joint probability distribution 
P n {x, k, e, /j,). This is the probability to end in position x with k reversals of direction 
after n steps. The random walk can be used as a simple model of a polymer with n 
atoms. Energy H is proportional to the number of kinks k and is given by H — —hk. 
This means that the polymer is totally folded at zero temperature. The position of 
the end point x measures the effect of an external force applied to the end point. 

In jS] we have calculated an exact expression for P n (x, k, e, fx). It is the product 
of two contributions 

P n (x,k,e,fi) = C n {x,k)F n (x,k,e,fi). (1) 

The function C n (x, k) counts the number of walks which start in the origin and end in 
X after n steps with k kinks. This function only depends on the final values of x and k. 
The probability that an n-step walk ends in x and has k kinks is denoted F n (x, k, e, //). 
This probability depends on the values of e, [i and the initial conditions. See figure ^ 
for an example of a walk with 6 steps and 4 kinks, ending in x = and starting to the 
right. There are two possible configurations Cg(0,4) = 2. The probability for both 
configurations is -^(0, 4, e, fx) = po(l — e) 2 (l — m) 2 ^- 
One has approximately 

P n (x,k,e,n) xip n (x,k,e,Li) = c n (x, k)f n (x, k, e, /i) 



f n (x,k,e,fi) = (1 -e)2(l -^e" ia + "^ n L° 
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The lattice parameter is denoted a. The approximate distribution p n (x,k,e, 
deviates from the exact distribution P n (x, k, e, /i) for two reasons. Firstly, P n (x, fc, e, ^) 
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has a small dependence on the initial conditions. Secondly, the expression P n (x, k, e, [i) 
is different for even or odd k. These effects are negligible when the number of atoms n 
is large. In the present paper, we only use the approximate distribution p n (x, k, e, /x). 
Define the following reparametrisation 

F=-Llni and = JL]n (1 " (5) 
Then the probability p n (x, k) (omitting the dependence on e and /z) can be written as 

p n (x,k) = 5&$-e-K- hk - F *\ (6) 
with the partition sum Z n approximately given by 

Z n ~ M" n/2 . (7) 

In this form, it is clear that p n (x,k) is a Gibbs distribution. We will show that the 
parameter (3 is the inverse temperature and F is the applied external force. 

The model depends on two parameters, the probabilities e and /j, (or the inverse 
temperature (3 and the applied external force F). Hence measurement of two 
observables is needed to estimate the values of all parameters of the model. In [5] 
we obtained exact expressions for the average number of kinks and the average end- 
position 

(k) = 2n K , (x)=na-— (8) 

These relations link the model parameters e and /i to the measurable quantities (x) 
and (k). 

4. Isotensial experiment 

The isotensial experiment is described in the grand canonical ensemble, with the end- 
to-end distance x and the energy H = —hk both fluctuating. The probability for a 
polymer with n + 1 atoms to end in x with k kinks is p n (x, k) and is given by ©■ The 
entropy is 

^S=- VV P „(^ *) In = -PHk) ~ (3F(x) + InZ n . (9) 

Hence, the free energy G equals 

G = J B-F(x)-T5 = -/3" 1 lnZ„. (10) 
From 1 = k p n {x, k) and © then follows 

° = ^I>^fc) = -£+^(z>- ^in^ ( n ) 

SC,fe 

Q = §fH^{x,k)=^{x)-^ f \uZ n . (12) 

Ur,fc 

With these relations one shows that the free energy G satisfies the relations 

SLpG = E-F(x) (13) 

^pG = - (x). (14) 
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Figure 2. Plot of F((x))a/h as a function of (x)/na. The value of 0h is 5 for 
the solid line and 2 for the dotted line. 



This proves that &b/3 is the inverse of the thermodynamic temperature T and that F is 
the applied force. Note that T should be positive. This gives the condition e + /j, < 1. 
By inverting (JSJ), one gets expressions for e and fi as a function of (fc) and (x) 

a(n - (fc)) + (x) 



a(n — (fc)) — (a:) 



na — (x) 

In combination with JSJ this gives 

= _J_ ln / (n- (fc))a+ (s) na-(x) 
2a/3 \ na + (a;) (n — (&))a — (x) 

In the latter expression (k) can be eliminated using 

(k) = - _ p X _ 2fflt fn - ^(x) 2 (1 - e- 2 / 3 ' 1 ) + n 2 a 2 e- 2 ^ 

The result is F as a function of (x) at constant temperature. It is plotted in figure 
for two values of (3h. 



(15) 
(16) 

(17) 

(18) 



5. Isometric experiment 

We proceed with the isometric experiment, which is described in the canonical 
ensemble, with fixed end-to-end distance x. Without loss of generality, we assume 
that x > 0. The probability distribution is the conditional distribution obtained from 
p n (x, k) by requiring that x has a given value 

_ p n (x,k) _ C n (x,k)e l3hk _ 1 k) J3hk /,QN 
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Figure 3. (A) Plot of the isometric force as a function of /3h, with n = 150 and 
x = na/2. An exact expression for this force is given by 1221 (solid line). After the 
continuum approximation we obtain 1231 (dotted line). (B) Plot of the adiabatic 
force as a function of j3h, with n = 150 and x = na/2. An exact expression for 
this force is given by IB. 11 (solid line). After the continuum approximation we 
obtain IB. 41 (dotted line). 



with z n = ^2 k c n (x, k)e^ hk the canonical partition function. The entropy is 

TT S = -5>»(fe|x)ln^4 = ~ Ph{k) +hXZn ' (20) 
/cb c n [x, k) 

The free energy is 

G = E-TS= -/5- 1 lnz„. (21) 

The mean force can be obtained by taking the derivative of G with respect to x. The 

latter is a discrete variable, so the force becomes 

, „. , . AG 1 A In Zn , , 

(F)(x) = — = --— 22 
Aa; p Ax 

We cannot calculate this discrete derivative in closed form. In the next section we 
proceed by assuming that k and x are continuous variables. Under this assumption an 
explicit expression for (F)(x) can be obtained. The differences between this expression 
and expression it22|) will be discussed in detail. See also Fig.|3J\. 

5.1. Continuum approximation 

Assuming that k and x are continuous variables (n, k and x large), one can replace 
summations by integrations. In |Appcndix A| we derive under this assumption following 
expression for the force 

WW =~i( r ,, dC i X,k) ) + l-Pn(n - x/a\x). (23) 



[3 \ Cn(x, k) dx I (3 a 
In the continuum approximation, the force consists of two contributions. The first 
contribution is due to the change of entropy, the second is a finite size effect (this 
term disappears in the limit n — > oo). 

Expressions 1221) and l|23|) differ in the limit (3 — ■> oo. Indeed, the limiting value 
of expression (|22|l is 

lim (F)(x) = - A lim (hk) = -h A (n - -) = -. (24) 
P— ►oo l\x /9^oo /Xx V a/ a 
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The limiting value of expression (|23|1 is 0. The force vanishes in the limit n — > 
oo because an infinite chain implies that the average number of kinks is infinite, 
independently of the end-to-end distance. This has the immediate consequence that 
for low temperatures the force equals zero. 

Mathematically, the origin of the difference between the force for a finite or an 
infinite chain, is the order of taking the limit n — > oo and performing the derivative 
with respect to x. This can be seen by taking the limit n — > oo of expression l)24fl 

- lim — lim (hk) = - ^ - — lim lim (hk) = 0. (25) 

We conclude that the order of taking limits influences the result for the force. To 
illustrate this we plotted expressions and as functions of (3 (with n and x 
large) in Figure^- The fi gure demonstrates that the continuum approximation holds 
for small values of f3, but not for large values of /3, because of the different limiting 
values. 



5.2. Adiabatic experiment 

If the experiment is carried out fast enough, the entropy is constant and does not 
contribute to the adiabatic force, which is given by (F) = —A(k)/Ax. The results 
obtained in this case are found in | Appendix B| Figure[3}3 shows expressions (JB.lfl and 
<|B.4fl as functions of (3. Again, the continuum approximation holds for small values 
of P, but not for large values of f3. Figure |3)3 shows also that the discrete force is an 
increasing function of (3. We can understand this by following reasoning. 

In the limit (3 — > oo, the isometric experiment and the adiabatic experiment 
coincide, because in this limit the free energy is equal to the energy. The average 
number of kinks equals the maximum that is allowed under the constraint of the fixed 
end-to-end distance x. Extending the chain with Ax, results in a decrease of the 
average number of kinks by Ax/a. So the force becomes h/a. For lower values of (3 
the average number of kinks decreases. Then an extension of the chain by Air can still 
be compensated by a decrease of the number of kinks, but also by a rearrangement of 
the kinks, without changing their number. So the decrease of the average number of 
kinks is smaller than Ax /a and the force is smaller than h/a. 

The difference between the isometric and the adiabatic force can most clearly 
be seen in the limit j3 — * 0. Figures and [3)3 demonstrate that the isometric 
force diverges in the high temperature limit, while the adiabatic force converges to a 
constant. This limiting behaviour can be calculated in closed form using the original, 
non-approximated expressions for c n (x,k) as defined in S|. The high temperature 
limit of the number of kinks is 

2" 



1- — 



Tl 

lim (k) = — 

The high temperature limit of the adiabatic force is 



(26) 



WF>(z) = --, (27) 
/3^o a na 

while the limiting value of the isometric force is 

, i , . 1 A , v-^ , ,n 1 , l + x/na 

WW - -0SkX>(*.*) - Ta^T-^a- (28) 

Stirling's approximation and the continuum approximation have been used to obtain 
the last expression. 
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x/na, <x>/na x/na > <x>'na 

Figure 4. (A) Plot of the grand canonical force (dotted line) as a function 
of average end-to-end distance and the canonical force (solid) line as a function 
of end-to-end distance (with n = 150 and j3h = 4). (B) Plot of the quotient 
(F)(x)/ f((x)) as a function of (average) end-to-end distance (with n = 150 and 
/3h = 4) — See tHJ. 



6. Comparison of the two force-extension relations 

Figure OK shows the exact force-extension relations for the grand canonical (|17J) and 
canonical ill'L'l) ensembles. As can be seen, they differ. In certain limits however the 
two force-extension relations do coincide. 



6.1. Thermodynamic limit of the force- extension relations 

Assume that (x)/na r* 1/2 and {k)/n m 1/2. Then one can approximate the grand 
canonical force by 

F({x))^f((x))-±-^. (29) 
ap na 

with f{x) defined by 

/(a;)= jLi n f ^-;fe)) a +^ ( 30) 



2a/? \ (n — {k))a — x / 

Figure0j3 shows the quotient (F)(x) / f{(x)). From this figure it is obvious that the 
canonical force and f((x}) are approximately the same for most values of (x)/na. 
The grand canonical force equals f((x)) up to a term proportional to (x)/na. So the 
force-extension relations of the grand canonical ensemble and the canonical ensemble 
are equivalent if (x)/na << f((x)). This inequality holds when (x)/na « 1/2 and 
(fc)/n«l/2. 



6.2. Physical explanation 

In the grand canonical ensemble, one can calculate the average value (x)(F) of the 
end-to-end distance for any given temperature T and external force F, using the 
probability distribution p n (x, k). In general, this average value does not coincide with 
the most probable value x mp of the end-to-end distance. In 0, Keller et al argue that 
the difference between the two experiments is a consequence of the difference between 
(x)(F) and £ mp . We will follow the lines of [S] to prove this statement for our model. 
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First consider the isotensial experiment. With a fixed force and a given 
temperature one can calculate (x)(F). After inverting this relation one obtains F((x)). 
The most probable value of the end-to-end distance can be calculated by solving the 
equation 



A 



0. 



(31) 



Based on numerical evidence, we assume that this equation has only one solution. 

Next consider the isometric experiment, with fixed end-to-end distance taken 
equal to x mp - Using formula ill' 21) one obtains for the average force 



k/h=0 



— x I a 



(32) 



In the continuum approximation following equation holds 

F((x)) = (F)(x mp ). (33) 

This is the same result as that obtained by Keller et al 6 . It means that the average 
force (F) measured in an isometric experiment equals to the fixed force F of the 
isotential experiment provided that the fixed extension of the isometric experiment 
is x mp instead of the average end-to-end distance of the isotential experiment (x). 
Therefore, the force-extension relations of the isometric and isotential experiment 



coincide only when (x) 
are negligible. 



This is the case in the limit of large n, when fluctuations 



7. Discussion 



In this paper we study two polymer stretching experiments. The theoretical 
description involves two different ensembles, one canonical, the other grand canonical. 
We have derived expressions for the force-extension relations in both ensembles using 
a simple model. 

The model is a one-dimensional random walk depending on two parameters. We 
showed in jSj that the probability distribution of the walk is approximately a Gibbs 
distribution, sec JBJ. In the present paper we ignore the small deviations from the 
Gibbs distribution. The model represents a polymer which has k kinks and end-to- 
end distance x. The averages (x) and (fc), or the corresponding control parameters 
external force F and temperature T, are used to determine the two model parameters. 

In the grand-canonical description the two observables k and x both fluctuate. 
The external force F is modelled as usual by adding a term —Fx to the Hamiltonian. 
The force F((x}) as function of the average end-to-end distance is obtained by inversion 
of (x) as a function of F. In the canonical ensemble, the end-to-end distance x is fixed. 
The probability distribution for this experiment is obtained from the grand canonical 
distribution by conditioning on x. The average force (F) is then obtained as minus 
the gradient of the free energy (or of the energy in the adiabatic case) . In this way 
one obtains the average force (F)(x) as function of the fixed end-to-end distance x. 
A peculiarity is that the result depends slightly on the order of taking the continuum 
limit and calculating the gradient. 

In general the two force-extension relations F((x)) and (F)(x) are different, which 
can clearly be seen in figure 0J\. In 6 , Keller et al argue that the difference between 
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the experiments is due to the difference between the average end-to-end distance 
calculated in the grand canonical ensemble and the most probable end-to-end distance 
of the grand canonical distribution. Our calculations confirm this point of view. But 
we find additional effects due to finiteness of the polymer chain. All these effects 
disappear in the limit of infinite long chains, so that the two experiments coincide in 
the thermodynamic limit. 

It is highly desirable to extend our work to two or three dimensions because the 
role of entropy is too limited in d — 1 . This is the reason why the present model does 
not exhibit a folding transition. 

Appendix A. Continuum limit 



Assume that x and k are continuous variables. Then one can write 



(F)(x) 



1 d 



In Zn 



1 1 dz n 
(3 z n dx ' 



f 

JO 



(A.1) 



[3dx 

The partition function becomes 

pn — x/a />oo 

dkc n {x,k)e pk = I dkc n {x,k)e /3hk Q(n-x/a-k) (A.2) 
io Jo 

with Q(x) — if x < and 6(2;) = 1 if x > 0. Calculating the derivative gives 

(FX*) = IJL ( I I" dkc n (x,k)e^ hk 5(n-x/a-k)- f"'^ dc "( x > k ) e phk\ 
P z n \a J J ox I 



1 dc n (x,k) 



1 1 < ,n 1 
= -3-Pn {n - x/a \x - - . 

pa p \ c n (x, k) ox 



(A.3) 



Appendix B. Adiabatic experiment 

If entropy is constant, force is obtained by taking the derivative of the energy 

AE 



(F)(x) = — = A(x) + B(x), 



with 

A(x) = 

hk 



En— x/a / 7 \ 
k =o c n (x,k)e 



fjhk 



En — x/a / j \ 
k=0 c n (x,K)e 



0hk 



' TT=o a ^ cn(x + a, k)e^ k J2lZo /a+1 ^{x - a, k)e^ k 

[c n (x, k)\ \ I 1 Ac n (x,k)\ I hk Ac n (x, k) \ I [c n (x, k)\ 



c n (x, k) j \c n (x, k) Ax 



B(x) = 7^Vn{n — x/a + l\x - a) 



c n (x,k) Ax 

n — x 
fc/h= 



Cn{x, fc) 



+ —p n (n - x/a\x + a) 



x ^ ^2k/h=o kc n( x + °j k)e@ hk 
J2Zo /a kcn(x-a,k)e^ k 



a Yr k ZZ' a+l c n {x-a,k)e^_ 



where we used the notation 



(B.l) 



(B.2) 



[cn(x, fc)J = - [c n (x + a,k) + c n (x - a, k)] . 



(B.3) 
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